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The study of the spin dynamics of a regular lattice of 
rare-earth ions in the paramagnetic phase of an insulator 
is not a subject of much interest in itself. However, ex- 
otic dynamics may has been seen at low temperature in 
geometrically frustrated magnetic rare-earth compounds, 
see Ref. [3| for an example and Ref . [2] for a general refer- 
ence. Therefore we find useful to establish the expected 
thermal behaviour of the crystal-field linewidth and the 
positive muon (/iSR) spin-lattice relaxation rate arising 
from the phonon-induced relaxation. 

A theory for the linewidth of crystal-field excitations 
in metallic rare-earth systems has been given a long time 
ago'^. The description of the single- ion dynamics was 
based on the memory function technique. Later on, the 
same problem was considered by Richards using a more 
conventional iterative method akin to the Rcdfields's, 
the derived system of equations being solved by Laplace 
transform^. 

Here we examine the contribution of the conduction 
electrons and phonons to the spin dynamics in the para- 
magnetic phase of compounds with a rare-earth sublat- 
tice. The crystal electric field acting on the magnetic ions 
will be dealt with. The interaction between the ions will 
be accounted for. 

Two experimental methods are commonly used to in- 
vestigate spin dynamics in magnetic materials: the neu- 
tron scattering and /iSR techniques. While the relation- 
ship between the measured intensity and the physical 
properties of the spin correlation functions is well doc- 
umented for the former, see for example Refs. [10] j it is 
not the case for the latter, at least for a lattice of rare- 
earth ions for which crystal-field effects are expected to 
be important. This is done in this paper, building on 
previous works dealing with the critical spin dynamics 
and magnon-induced relaxation as probed by /iSR, see 
Ref. and references therein. 

The organization of this paper is as follows. In 
Sec. m we investigate some basic properties of the spin- 
correlation tensor. In Sec. |TT] we write the /xSR spin- 
lattice relaxation rate, A^, in terms of spin-correlation 



functions. In Sec. lIIII these functions are described for the 
case of a single rare-earth ion, following the method of 
Richards'*. In the following section. Sec. IIVI we give the 
expressions for the conduction electrons and phonons cor- 
relation functions. The theory develops in the previous 
sections is applied for simple crystal-field level schemes 
in Sec. |Vl The interaction between the magnetic ions 
is taken into account in Sec. |Vll We start by using the 
random-phase approximation (RPA). Since this approxi- 
mation is found to breakdown for a limit of experimental 
interest, a phenomelogical modification is proposed. 



PROPERTIES OF THE SPIN-CORRELATION 
TENSOR 



In this section, we study some basic properties, 
mostly resulting from the point symmetry at the 
rare-earth site, of the symmetrised spin-correlation 
function Aq^(cj) = {{Ja{t)Jp}) exp{iujt)dt, where 

2 {{Ja{t)Ji3}) = {Ja{t)Jp) + {JpJait))- Ja is the a Com- 
ponent of the total angular momentum of a magnetic ion 
in the compound of interest. We explicitly assume J to 
be a good quantum number. This is true for a rare-earth 
ion. 

Our first task is to examine the structure of A(aj), the 
tensor elements of which are the functions Aap{Lo). This 
tensor is written in the reference frame {x, y, z} attached 
to the crystal lattice. In fact, because the crystal electric 
field potential is expressed in terms of Stevens operators 
which are polynomials of J+ — Jx -\- iJy, J- = Jx ~ iJy 
and Jz, we shall work with the reference frame {+, — , z}. 
It is only in this representation that the Stevens operators 
display nicely their symmetry properties. 

In this study we shall split the total Hamiltonian of the 
system under study into two parts. We write Ti. = Ti.o + 
Til where Tii describes the perturbation of the system 
relative to Tio- In addition, we take Tip as the sum of two 
commutating operators. Tip = '^cf + Wl, where the first 
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term is modelling the crystal-field energy levels of the 
rare-earth ions and the second the lattice. We restrict 
ourselves to the analysis of zero-field measurements. We 
shall assume in this paper [Hl, Ja{t)] = 0. 

Physically the relaxation of a spin correlation function 
of a single rare-earth ion is driven by the crystal-field 
transitions induced by the coupling of the crystal field 
with the lattice. Neglecting this coupling, i.e. assuming 
Til = 0, we determine which spin correlations are com- 
patible with the point symmetry at the rare-earth site. 
If a correlation vanishes when the coupling is neglected, 
it cannot be large even if it is taken into account. Hence, 
we study the symmetry properties of (Ja{t)Jp)^, where 

— Tr {poexp {iHot/h) JaCxp {—iHot/h) J/3} . (1) 

We have introduced the density operator po — 
exp (— /3Wo) /Zq where Zq is the partition function of 
T-Co, P ^ l/(fcBT) and Tr{A} stands for the trace of A 
over the quantum states of Tio- Since [Hl, Ja(t)] = 0, 
[Hl, Ja] ~ and thus 

= Tr{pexp{iHcFt/h)Jaexp{~mcFt/h)Jp}. (2) 

p is the density operator for the crystal-field levels. We 
use the simplified notations p and (• • •) rather than pcF 
and (• ■ ■)fjp because there is no risk of confusion. 

We are going to show that the only tensor elements 
which may not be zero in most physical cases are A^^i'-^), 

A_| {(jj) and A ^{uj). To proof, for example, that A^z{^) 

is usually incompatible with the point symmetry at the 
rare-earth site, we need to show {J^{t)Jz) = 0. We have 

{J-{t)J.) 

= Tr {pexp {iHcFt/h) J_ exp {—iTicFt/K) Jz} 

= Tr {pexp {iHcFt/h)IJ^Iexp {-iUcFt/h) (3) 

We have inserted the unit tensor X twice. Hcf is invari- 
ant under the operations of the point symmetry group 
at the ion site. Let us assume that one of the symmetry 
operation of the group is a rotation of angle g = 27r/n 
about the z axis, i.e. the z axis is a n-fold symmetry axis. 
Mathematically this invariance means 

[HcF,7^,(g)] =0, (4) 

where TZz{q) = exp(—igJz) stands for the rotation oper- 
ator about the z axis of angle g; see for example Ref. Q. 
Recalling the relation 

TZzig)J-TZ^^{g) = exp{-ig)J_, (5) 

and Eq. |4l starting from Eq. [3] we derive 

(J_(t)J,) = exp(-ie)(J_(i)J.). (6) 

We have used X = TZ~^{g)TZz{g). The only solution of 
this equation with g = 2n/n and n > 2 is {J^(t)Jz) = 0. 



In the same way we find {JzJ-{t)) = 0. Our conclu- 
sion is that the existence of the symmetrised correlation 
function ({ J_(i) J^}) is inconsistent with the z axis being 
two-fold or of higher symmetry. Using the same proce- 
dure we find five other elements of A{uj) to be zero when 
n >3, because ({J_(f)J_}) = only if n > 3. However, 
this tensor is not diagonal in {+, — , z} since the three non 
vanishing functions are A.| (w), A ^{lo) and A^zi^o). 

Before diagonalizing A(a;), we first establish some 
general properties of the symmetrised correlation func- 
tions. As {Jait)Jp)* = (4(-i)4) and {Jc^Jpit))* = 

{JpJa{~t)), from the definition of Azz i^^), A^ (uj) and 

A ^ (uj) we infer these functions to be real. In our nota- 
tions At is the adjoint of operator A. Since in addition 
{Mt)Jp)* = (44W) and {Ja^Jpit))* = {4{t)Jl), we 
derive A_| (w) = A |_ (—1^). Therefore the present dis- 
cussion leads to restrict the number of correlation func- 
tions to be computed to two: Azz i^) and, e.g. A_| (uj). 

A{uj) is diagonal in a frame denotes as {+',—', z}. We 
shall need to specify explicitly only the eigenvalues. They 
are A+/(w) = A+/+/(i:j) = ^A+_(cj)A+_(-a;), A_,{lj) = 
A_/_/(ti)) — — A_(-/(w) and Az{lli) — Azz{(^)- One easily 
shows that A^ (w)A^ {~'^) ^ 0- 

Let us introduce the auxilary function 

/oo 
{Ja{t)Jp)exp{iujt)dt. (7) 
-00 

The condition of detailed balance gives the relationship^j^ 

flapiuj) = exp {(ihuj) Q.pa{—bj). (8) 

This implies 

A+_(w) = [l + exp(-/3Sw)]f7+_ (w)/2, 
Azz{^) = [l + exp(-/3M]f^..H/2. (9) 

Hence, once VLzz (w) and fi^ {lo) are computed, A^ [lo) 

and Azz (w) can be evaluated. 

Note that the inelastic neutron scattered intensity is 
directly proportional to VLap{uj)\ see for example Refs. 
Q. The so-called scattering function Sap{i-o) is within a 
constant of proportionality equal to rict^(tj). 

II. ^iSR SPIN-LATTICE RELAXATION RATE 
AND SPIN-CORRELATION FUNCTIONS 

The /iSR techniques are based on the detection of the 
positron produced during the muon decay. Of particular 
interest is the longitudinal field technique for which the 
initial muon beam polarization, the external field Bcxt if 
any, and the axis of the positron detectors are all parallel 
to the same axis denoted as Z. Therefore that axis is 
of special importance. We define it by its polar angles 9 
and f in the reference frame {x, y, z} introduced in the 
previous section. We note that the /iSR techniques probe 
only the magnetic field at the muon site and therefore 
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the magnetic field correlation tensor $ (w^) at an energy 
hbjf^. uj^ is proportional to the local field at the muon site. 
Since the difference between that field and Bext is usually 
very small^, we write cj^ — 7^i?cxt where 7^ is the muon 
gyromagnetic ratio (7^ = 851.6 Mrad s~^ T~^). (t^^) 
is also defined in {+', — ',z'}. 

A longitudinal field /iSR experiment consists of mea- 
suring the relaxation function which is, to a good ap- 
proximation, an exponential function in the case of in- 
terest here, e.g. for paramagnets with a regular lattice of 
rare-earth ions, and therefore is characterized by a single 
parameter, the spin-lattice relaxation rate Xz- Implicit 
is the assumption that the motional narrowing limit ap- 
plies, i.e. the fluctuations of the rare-earth moments are 
rapid. Az is directly related to $ (w^): 

= YE^7M^>v')'i>/37K), (10) 

where now {7,/?} — {+',—',2}, the matrix elements 
Ljfj{6,ip) can be com put ed from the material given in 
Appendix C of Ref. \ld\ and (ujf^) is an element 
of $(0;^). We have introduced the Fourier transform 
^l^j {ui) = ^fjj {t) exp{iujt)dt, a definition differ- 

ent from the one used previously^i. The total angular 
momenta of the magnetic ions interact with the muon 
spin through the dipole interaction and, in addition, the 
hyperfine interaction which accounts for the transfered 
magnetic field for an insulator and the conduction elec- 
trons for a metallic compound. This induces a magnetic 
field at the muon site with correlation functions such that 

a i,i 

(11) 

i and j run over the magnetic ions. Ka-ij {uj) is a spin- 
correlation function between sites i and j. is the 
coupling tensor labeled by the vector which points to 
site i and the origin of which is the muon site. Explicitly 
we have 

= -{^^-^-^)+Hr.S^''- (12) 

Hr- is the hyperfine constant for the coupling to the ion 
i. For simplicity we assume all the magnetic ions to be 
equivalent: they sit on a Bravais lattice, no is the mag- 
netic permeability of free space, g the spectroscopic fac- 
tor, /iB the electronic Bohr magneton and v the volume 
per magnetic ion. 

Two approaches are now possible. If the interaction 
between the total angular momenta of the magnetic ions 
is negligible, we can restrict ourselves to the behaviour 
of a single momentum. Otherwise their interactions has 
to be taken into account. 

We first consider the case for which Aa-^ij (w) = 
Aq {lu) Si,j. The intersite correlations are therefore negli- 



gible and thus 

^'^-(^^ - 0'^EA"(-)EGi?rG-. 

a i 

(13) 

Hence, based on our discussion in Sec. [H Xz is found 
to be the weighted sum of two spin-correlation func- 
tions for a compound with at least a three-fold symme- 
try axis: they are A^ (w) = K^z (w) and, e.g. A_|_/(w) = 
•\/A-| (i-ij)A^ {~^)- Obviously, additional symmetry el- 
ements of the point group at the rare-earth site may re- 
duce the number of correlations to be computed to only 
one. 

It is convenient to use the space-Fourier transform 
Ka{(\,oj) if the intersite correlations are appreciable. Re- 
ferring to Ref. , 

a. 

(14) 

with the definitions ^''^"(q) = G''"(q)G"''(-q) and 
G/5«(q) = G^^ exp (iq • r,) . V is the volume of the 
sample, i.e. V — Nv where A'' is the number of magnetic 
ions in the sample. The integration is over the Brillouin 
zone. Again, two spin-correlation functions are sufficient 
to model Xz for a compound with at least a three-fold 
symmetry axis. As pointed out above, the number may 
be reduce to only one due to the presence of additional 
symmetries at the rare-earth site. 

III. SINGLE-ION CORRELATION FUNCTIONS 

Q^p (lo) in the reference frame {+, -,z} 

In this section we give formulae for Qzz (w) and 

fl^ (a;). The interaction between the rare-earth ions is 

neglected. The correlation functions are evaluated with 
the Richards's method^. In fact, it is the Laplace trans- 
form of Ja{t) which is considered there. P. M. Richards 
has only presented the result of his work, focusing on the 
electronic relaxation. Here we sketch the derivation and 
apply the result for two mechanisms which can drive the 
crystal-field transitions: the conduction electrons and the 
strain field. 

The interaction representation is used. This requires 
to introduce the two new operators, Ja{t) and HKt) (see 
for example Ref [3]), such that 

J*{t) = exp{—iHot/h)Ja{t)exp{iHot/h), 

n{{t) = exp{-inQt/h)Hicxp{inQt/h) . (15) 

Note that J*{t = 0) = Ja(t = 0). 

Using the Heisenberg evolution equation for Ja{t), 



exp {iHot/h) °^ ' exp (-iHot/h) 



with 



dt 



dt 



(16) 



(17) 



Solving this last equation for J*(i), 



We are basically interested by the matrix elements of 
Eg. 1161 For their evaluation, we introduce the expression 
for J* (t) in the commutator of Eq. \17\ dropping the first 
term on the right hand side which is time independent. 
We note the kets of TicF as {\m), \n), ■ ■ ■} with, for ex- 
ample, HcF\m) = Hu!m\m). Denoting UJrnn = - ^^n 



and recalling the assumption [Hh, Ja{t)] = 0, we derive 
dJait) 



dt 



- iuj.mn {m \Ja{t)\ u) 



I + 11 + III + IV, 



(19) 



where the four terms labelled by IV have similar 

mathematical structures. 

We need to specify Hi at this junction. We focus first 
on the crystal-field fluctuations induced by conduction 
electrons. The exchange interaction between a rare-earth 
ion and the conduction electron can be written^ 



Hi = Hel = -2/ex(.9 - 1)J • S, 



(20) 



where /ex is an exchange integral, g the Lande factor and 
s the electron spin operator. Then, for example. 



fi2 



Hi exp 



HiJa{t') exp 



(t-i) 



dt' 



4/cl(g - 1) 

P:9 



2 /.t 

- / y^X!*'^ ^™ 1-^71 P) exp [iLOpn (t - t')] {t - t') (p I J^' I (?) {q \Ja{t')\n) dt', 
"'O ,„ „ „ , 



p,<i in 



exp [iujpn {t - t')] /oi [- {t - t')] (q \Ja{t')\ n) dt', 



(21) 



with G|i,,„„ = [4llig - l)yh^] (fc \J,\i) {m \J^\ n). 
In the last step we have performed an average over the 
electronic states assuming the electronic correlation to be 
isotropic and given by fe\{t) = {sx{t)sx)^y 

Before discussing the other terms for the electronic 
relaxation, we pay attention to the / expression for 
the phonon-induced relaxation. From numerous electron 
paramagnetic resonance (EPR) studies, see for example 
Ref. |lj|, the phonon-induced relaxation between elec- 
tronic energy levels is known to be driven by the modu- 
lation of the crystal field. We expand that field in powers 
of the strain and limit the expansion to the linear term 
denoted as Hsp ■ It accounts for the spin-phonon interac- 
tions, i.e. the magneto-elastic coupling, and represents 
an additional electric potential generated by the lattice 
vibrations. To determine its form, we should study in de- 
tails the vibration-induced distorsion of the surrounding 
of the ion of intereslji^. Here we follow a much simpler 
path. We take an average strain e and ignore any direc- 



tional properties. This results to 

Hi = Hsp - eV (22) 

To second order in J, or more practically using the 
Stevens operators (these operators are listed in dif- 
ferent references, e.g. Refs. (l3fl6j ). we have 

V = goOl+giOl+g.iO^^ +g20l+g.20^\ 

(23) 

where the five gn are unknown real parameters. Only 
the operators in this sum which break the point sym- 
metry at the rare-earth site are effective in inducing the 
crystal- field transitions. Practically, to limit the number 
of parameters, one may have to select one or two out of 
the sum in a somewhat arbitrary fashion. Focusing again 
on the / term. 



iujpn {t - t')] /ph [- (t - t')] {q I J„(t')| n) dt', 



(24) 
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with Gf,^^^ = {k\V\i){m\V\n)/h'' and /phW = 
(e(<)e)pjj. Here the average is done over the phonon bath. 

It is obvious that the / term for the electronic and 
phonon relaxation processes have the same mathemati- 
cal structure. This is also true for the three other terms, 
that is, //, /// and IV . Therefore we need to solve Eg.fTOl 
for Ja (t) with the four terms on the right hand side hav- 

I 



We denote g{uj) = f{t) exp(icji)dt and h{uj) — 
ir exp(iwi)dt. Gki^rnn and f{t) stand for Gf^ „^^ 

and fci{t) or G^^ and fph{t), depending on the relax- 

I 



As argued in Sec. llVi it is a good approximation to ne- 
glect the imaginary part of ft, (a;) and g(uj). We note that 

h{uj) ^ g{-u}) ^ exp{-l3huj) g{Lj). (27) 

The inclusion of the term {Ja)5mn guarantees that 

I 



Sl+_(cj) = 



c.c. stands for the complex conjugate of the expression 
on the left of this symbol. Zqf — J2m {~l3Em) is the 
partition function for the crystal-field level scheme of a 
rare-earth ion and the Za,mn functions are obtained by 
solving the system of coupled (2J -f 1)^ linear equations 
given at Eq. [26l 



ing the mathematical structure of Eq. [21] (equivalent to 
Eq. [M)) . Since each term is expressed as a sum of con- 
volution products, a Laplace transform technique is well 
suited. That TLi factorises into lattice and crystal-field 
potentials is seen to be a key ingredient which allows to 
use the Laplace transform technique. With the notation 
Za,mn{s) = /(,°° exp(-st) (to| Jq (t) | n)dt , 



(25) 



ation mechanism under study. Setting s — e + iu) with 
e ^ 0+ in the previous system of equations. 



(26) 



{Ja{t)) relax toward its thermal equilibrium value. This 
is further discussed for the Rcdficld theory by, for exam- 
ple, Shchteri^. 

In fact, as explained in SecU Q.^ (w) and Q.^^i'^) are 

required. From their definitions. 



(28) 



IV. CONDUCTION ELECTRON AND PHONON 
CORRELATION FUNCTIONS 



In this section we evaluate the correlation functions for 
the conduction electrons and phonons. 

Turning our attention to the electronic correlation, we 



+ ^ Gmp,qn [h {is + UJpn) + g {is + UJmq)] Za,pq{s) + (m \Ja\n) . 



i {u! — LUmn) Za.mn{i^) + Gmp,pq h {—UJ + LOpn) Za,qn{i^) + Gpq^qn g {~U! + UJ,nq) Za.rnp{i^) 
- ^ Gmp,qn [h {-UJ + UJpn) + {-^J + UJ^iq)] Za,pq{iuj) = (m \Ja\n) - {Ja)5mn- 



P,<1 



—^—V'exp (-/?£;,„) Z+,™„(-iw) (n|J_|rn) -|-c.c., 

ZCF 

m.n 

-J—y^ exp{~PE,n) Z^,„in(-iuj) {n\J^\m) + c.c. 
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need to compute 

/>oo 

5ci('^)= / /ciWexpM)dt. (29) 

JQ 

We shall denote cj^^ and Ck,a the creation and annihila- 
tion electron operators of wavevector k and spin a and 
He the Hamiltonian of the electron bath. The electrons 
are assumed to be free-like. We use the fermion repre- 
sentation of Sx, i.e. 

= ^ ('^^M^k.T + cJ'jCka) • (30) 

k',k 

I 



Without lost of generality, the ion is taken to sit at the 
origin of the coordinates. /ci(i) reduces to the sum of 
two quantities which occur to be equal. Let us focus on 
one of them: 



{^CXp [iUet/h) 4,_|Cka exp {-iUet/h) cJ^^CqJ^^^ 

k,k'' q.q' 

= I] [* (^k',T - ^k,i) ^] (ctjCk,iC+ _^CqJ^^^ 

kjk' q.q' 

" ^ZZ^^P[*(^k',T - Wka)^]'^k'.qV,k(ck',|Cq,T)^j(ck,iCq',i)^l ^ ^1 " "^k) t] riq (1 - rifc) . (31) 

k,k' q,q' k,q 



riq is the Fermi distribution function. The second line 
results after inserting in the first line the identity oper- 
ator between the first two fermionic operators and using 
the Heisenberg equation. An Hartree-Fock decoupling is 
done and the spin dependence of the energy is neglected 
at the third line. The imaginary part of gc\{'^) leads to a 
slight energy shift and the real part describes the crystal- 
field relaxation. We shall neglect the shift and thus we 
only keep the real part of gc\{uj). It is probably difficult 
to detect experimentally the energy shift. We compute 

5ei(u;) = ^h[N{E^)f- ^ (32) 

2 1 — exp (—fjnoj) 

N (E-p) is the density at the Fermi level per spin. 

We now evaluate the phonon correlation function. It 
is convenient to express e in terms of its Fourier compo- 
nents. Since the rare-earth ion can be taken at the origin 
of the coordinates (see for example Ref. ^17|). 




2AfMwk 



(«k - 4) 



(33) 



N is the number of rare-earth ions in the compound and 
NM its mass, and ak are the creation and annihila- 
tion phonon operators. Note that only a single branch of 
phonons is accounted for. It is a simple matter to derive 



/ph 



2MuJk 



6k is the Bose distribution function. Because of energy 
conservation during the scattering process, the second 
term cannot contribute to the relaxation when a; > 0. 
We compute 



5ph(w) 



IV ^ 2Mcjk 

k 

3 1 



(6k + l)S{u - Wk) 



47r Qv^h? 1 — exp {—(3huj) 



(35) 



X [exp (-iwki) (6k + 1) + exp {iuj-kt) 6k] . (34) 



The imaginary part of g-p\^{ijj) is again neglected. The fac- 
tor 3 is added to account for the three phonons branches. 
The k sum, which extends over the Brillouin zone, is per- 
formed assuming lo\s^ = Vsk. g is the crystal density and 
Vg the sound velocity. 



V. MODEL COMPUTATIONS OF THE 
CORRELATION FUNCTIONS AND THE ^iSR 
SPIN-LATTICE RELAXATION RATE 

In this section we present analytical results obtained 
for simple rare-earth level schemes, our purpose being to 
understand the implications, on the correlation functions 
and /iSR relaxation rate, of the crystal-field-induced re- 
laxation by electrons and phonons. Here we shall focus 
our analysis on Q,zz {^)- As mentionned in Sec.[Tll is 
the weighted sum of two spin-correlation functions for a 
compound with at least a three- fold symmetry axis. We 
shall assume, for simplicity, Xz to be proportional only 
to A^^ (w = 0) = (w = 0). 
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We first discussed a result obtained from the itera- 
tive method detailed here against a published one derived 
from the memory function method that was applied to 
the electronic relaxation of a Ce'^"'" in cubic symmetry^. 

The six crystal-field levels of Ce^+ split into a dou- 
blet T^ and a quartet Fg located at such a high en- 
ergy than it can be neglected for our purpose. We 
label the two states of T-^ as |T) and |2}. They are 
|I) = (l/6)i/2| 5/2 } - (5/6)i/2| - 3/2 ) and |2) = 
(l/6)i/2| -5/2}-(5/6)i/2|3/2). We compute (1|J.|2) = 
(2| J^|l} = and (2| J;,|2) ^ J^|l). It is easily shown 
that ^2^12 = .^2,21 = and Z^ai = -Zz,22- Therefore 
we only need to study Eq. [55] for m = n = 1 . We expect 
rizz (to) to be a Lorentzian function. This is obtained 
setting uj — in the arguments of the h and g functions. 
This low-energy approximation will be always done in 
this section for simplicity. Obviously, it can be easily 
overcome. We derive: 

luZz^ii + AKiiuj = 0)Gf2^2iZzM - (36) 

The relaxation arises from a Gke.mn term that satisfies, 
as expected, the secular condition ojm — ^mrr^- We 
compute 



9r, 



z,cl 



2 u;2 



with the Korringa la^ 



M. 



kBT 



(37) 



(38) 



grj = 2(2|Jz|2) = 5/3 is the spectroscopic factor of T^. 
This result was published by Becker and collaborators^. 
From a straightforward computation. 



A..(w = 0) 



1 



8vr [I^^{g - 1)N{Ef)Y ksT 



(39) 



This model predicts Xz oc 1/T. Such a behaviour was 
first reported for ErAl2^. According to Eq. [311 the 
linewidth is proportional to hei{i-o — 0)- If we had consid- 
ered the phonon relaxation mechanism, we would have 
found the linewidth to be proportional to hph{i-o — 0). 
Hence, F^ ph = 0, i. e. in contrast to the electronic 
case the coupling of a doublet to the strain field does not 
broaden the quasi-static scattering. 

The work reported in this paper has been originally 
motivated by the observation for some geometrically frus- 
trated magnetic materials of a strong temperature de- 
pendence of Xz at high temperature and sometimes a 
plateau at low temperature for measurements performed 
in zero field^'^"'^^'^^. These materials are insulators, and 
therefore only the phonon mechanism can play a role. 

The physics at high temperature can be understood 
with a simple crystal-field scheme of three levels at en- 
ergies (2|Hcf|2) = (5, (3|Hcf|3) = A with A > 5 and 
zero. The ground state is the |1) state with energy equal 



to zero. Since our interest is on the low-energy physics, 
we study Eq. [^D with m — I and n ^ 2. In addition, 
assuming 5 negligible, 



(*^ + riS)^-.i2 = (i|J.|2), 



where 



.(12) 



V(A/n) (CP^^gi 



"^23,32 



(40) 



(41) 



The contributions to this equation come only from the 
two terms that satisfy the secular condition ujki = tUmn 
for G5^£^„. The other terms oscillate with an angular 

frequency lua = A/^ ~ lO^^j-ads"^ for A/fce 100 K, 
and therefore average out to zero in the time scale of 
the measurements. This simple model predicts a quasi- 
elastic scattering. Assuming further f3A ^ 1 gives 

f 12) 

F^ ph ^ exp(— /3A) and therefore an activated behaviour 
for Xz{T) is found. Components such as Zz,3i are not rel- 
evant for the quasi-elastic response and Xz since they are 
strongly inelastic. This argument is valid if the tempera- 
ture is not too large so that Z^^si is a Lorentzian function 
of UJ and so the contribution of Z^^si at w = is negli- 
gible. The relaxation we are discussing is Orbach-lik&i^. 
It involves a real two-phonon scattering: one phonon in- 
duces a transition between jl) and |3) and the other be- 
tween 1 3) and |2). It has been identified experimentally 
in /iSR data^S.. This scattering occurs only if phonons of 
the proper energy are available. For a compound with 
a Debye phonon spectrum, this means that the excita- 
tion energy A cannot be larger than fceOo where Qd 
is the Debye temperature. To complete our analysis 

(12) 

of the high temperature regime, we note that F^ is 
found to be proportional to T when /3A ^ 1. Therefore 
Xz (X 1/T when the thermal energy exceeds the crystal- 
field energy. Amusingly, the Korringa relaxation gives an 
identical temperature dependence. 

The same crystal-field scheme can be used to study 
the physics at low temperature. The high-energy level is 
now neglected. Focusing again on the low-energy sector, 
Zz ii and Zz^22 are the relevant components because S is 
no more negligible. They obey a system of two coupled 
linear equations: 

iujZz.ll + 2^ph(^/?l)G^2 21 (^2,11 - ■2'2,22) 
= {l\Jz\l}-{Jz}, 

'ii^Zz.22 ~ 25ph(5/?l)Gi2,21 iZz,ll - Zz,22) 

= {2\Jz\2)- {Jz). ' (42) 

This system can be solved for Zz^ii and 2^^.22 and Qzz{^) 
evaluated. Here we simply consider two limits. If /3(5 ^ 1, 
Zz,ii controls the quasi-elastic linewidth. Since hp^ and 
{{l\Jz\l)~ (Jz)) vanish, the width of the quasi-elastic line 
vanishes also. Hence, the coupling of the strain field to 
the crystal field does not induce a linewidth. In the limit 
/3(5 ^ 1 , Zz^ii and ^3^22 have an equal weight in the 
expression for Qzz{^)- Since they are of opposite sign. 
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f^zz('^) = 0. Hence, in the high-temperature hmit of 
the two-level system the magneto-elastic coupling does 
not induce a linewidth. However, as our study of the 
three-level system has shown, the energy levels located 
at energy above S may contribute to the relaxation. 

In summary, in this section we have found that the 
magneto-elastic coupling can effectively explain the acti- 
vated behaviour of Xz found at relatively high temper- 
ature for some insulators. The same coupling gives a 
vanishing linewidth at low temperature. However, other 
mechanisms can contribute to it. Such a mechanism is 
unveiled in the next section. 



VI. SPIN-CORRELATION FOR A REGULAR 
LATTICE OF PARAMAGNETIC IONS 

Here we study the effect of the magnetic interactions 
between the ions located in a regular crystal lattice. 
Hence, rather than Hq = Ti-cF + with T^l = ^ci 
or Hh = "Hph , we write Ho — Hcf + "Hn + Hh where 
Hu describes the ion-ion magnetic interaction. We shall 
limit ourselves to the two-ion Heisenberg Hamiltonian, 
i.e. 

= - J(*j)J.-Jj, (43) 

where J^(ij) = when i — j. Introducing the following 
relation in the previous equation, 

q 

we derive 

q 

We have used the relation 

Jq^^-iicxp{~ic[- i) . (46) 

i 

We recall that the magnetic ions are taken to sit on a 
Bravais lattice for simplicity. As in Sec. [TTl it is usually 
convenient to work with continuous rather than discrete 
wavevectors. This leads to 

nn = -vj J(q)J(q) • J(-q) (47) 

with 

J(q) = ^ J(ij)expHq-(i-j)]- (48) 

i 

From now on we shall work with the reference frame 
{+',—', z} where the spin-correlation tensors are diag- 
onal. To simplify the notation and because {to) is 



simply (w) (see Sec. HI]), we shall only consider that 
correlation. Rather than Az{w) we have to evaluate the 
wavevector dependent symmetrised correlation function 

/oo 
{{Jz{q, t)Jz{-(i, 0)}) exp{iujt)dt, 
-OO 

(49) 

where 2({J,(q,i)J,(-q,0)}) = (J,(q,i)J,(-q,0)) + 
(J,(-q,0)J,(q,t)). 

In general the available theories does not consider the 
correlation tensor but rather the generalized susceptibil- 
ity tensor. The fluctuation-dissipation theorem provides 
a relation between the elements of the two tensors, see 
for example Ref. [1^. Denoting the imaginary part of 
Xz(q,w) as XTO{xz(q,w)}, 

A^(q,w) = 3 coth (^^] Jm{xz{q,i^)}- (50) 

Therefore we need to compute Xmlxzil,^)}- 

A practical theoretical method to describe the effect of 
the intersite correlations on the susceptibility is the RPA 
which leads to the following simple formula^'! 

X(q,c.)-|x-^^^^x(^)j 'x(^). (51) 

{A}^^ stands for the inverse of tensor A. xi'^) is the 
susceptibility tensor for a magnetic ion in the lattice, the 
intersite interactions being switched off. The elements of 
the tensors xi^) ^^'^ ■^i'^) related. Using again the 
fluctuation-dissipation theorem, 

I-{x.(-)} = ^^tanh A.(.). (52) 

The real and imaginary parts of Xz(^) are linked through 
the Kramers-Kronig theorem: 

neixzic.)} = -'-V r ?!!!M!^d.. (53) 

J-oo 

V stands for the principal part of the integral. Thus, 
once tensor elements such as t^z{^^ £^re known, the ele- 
ments of xkiS) can be computed using Eqs. [52] and [53l 
Then, X■m\xz(s^.^ i^)} can be determined from Eg. 1511 Fi- 
nally A2(q, w) is obtained recalling Eq. [50] Within our 
conventions x(q, ^) and xf\Lo) are diniensionless. 

Since the tensors are diagonal in the reference frame 
{+', — ', z}, Eq.[ni]can be solved easily for x(q, Com- 
bined with Eq. Ho] 
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A^(q,tj) = 



hv 



coth 



Xm 



1 - 



2J(q) 



-PC 



A.(«) 

OO UJ — U 



tanh 



(^) 



2.7(q) 



tanh(^) A,(co) 



2 • 



(54) 



Therefore, \i Kz{lo) and J'(q) are known, it is possible 
to compute numerically (q, uj) and \z for a given muon 
site using the material of SecHH We can proceed further 
analytically for the elastic case, i.e. when a; = 0. For 
definitiveness, we assume h.z{Lo) to be Lorenzian, i.e. 



w2 + r? 



(55) 



We derive 



A,(q,c^-0) 



= 0) 



4(Jg).7(q) 1 TOO 



if- 1^ tanh ( ) dx 



c{a;2 + l) 



1 r 

TT J — 



■ tanh f ) da; 



(56) 



To understand this expression we consider the conse- 
quences of a magnetic phase transition. That transition 
is characterized, among other properties, by its critical 
temperature Tc and a wavevector qc which sets the rep- 
etition length of the magnetic structure for T < Tc- At 
Tc the generalized susceptibility diverges for q = qc, and 
so does A^(qc,w = 0). From the previous equation 



h _ 4(jf ) 1 

Fz TT 



J(qc) 



1 



■ tanh 



dx. 



(57) 



Let us study I{a) = x{x^+i) tanh(aa:)da;. We eval- 
uate numerically /(I) = 2.08. This is clearly smaller 
than r°° -2^., dx = n. However, smaller is a better is 
the approximation of linearising the function tanh(aa;). 
In our case a = {HT z) / {2kBTc) . Physically, we expect 
HTz — fceT'c for a conventional phase transition. Assum- 
ing the linearisation to be valid, we get 



J(qc) = 



(58) 



This is consistent with the mean- field result. 

We known that Xz in zero field is given as a sum over 
the Brillouin zone of Az{q_,uj = 0). Approaching the 
critical point, the sum is expected to be dominated by 
Az{q^c,^ = 0) since that correlation diverges at Tc. How 



does it diverge? According to our RPA result of Eq. [56l 
upon linearising the tanh function, we derive 



Az{qc,u; = 0) = A^(cj = 0) 



T 



T-Tc 



(59) 



While Eq. [55] is obviously not valid in the critical regime 
since the RPA is a kind of mean-field approximation, it 
should provide a good description deep in the paramag- 
netic regime. Note that Az{q^c,^ = 0) reaches a finite 
temperature independent value at high temperature. 

Now we study Az{q,uj = 0) in the limit T :$> Tc 
{Tc ~ Tc for a ferromagnet and Tc = Tn for an anti- 
ferromagnet). We find: 

A,(q,c. = 0) = A,(c.^0)(^l + 2|Ml£^. (60) 

Because J^{qc) > 0, the sign of the thermal correction 
at high temperature is given by the sign of J'(q). Re- 
ferring to Eq. 271 positive and negative J{q) values cor- 
respond to ferromagnetic and antiferromagnetic interac- 
tions at wavevector q, respectively. This means that 
a ferromagnetic (antiferromagnetic) interaction increases 
(decreases) the correlation function Az{q,uj — 0). This 
result was derived a long time ago, but the wavevector 
dependence of the interaction was overlookedi^i^. The 
1 /T dependence of Az (q, w = 0) is a typical pair correla- 
tion effect at high temperature. 

The result of the RPA approximation is quite sim- 
ple, in particular for Az{q,uj = 0). When Az{u! = 0) 
is determined, A2(q, uj = 0) can be evaluated by dividing 
Az{uj = 0) by a proper denominator. However, in Sec. fVl 
we have found A^ (w = 0) to vanish if the thermal energy 
is much larger than the crystal-field energy Ecf when 
the magneto-elastic coupling drives the crystal-field re- 
laxation. In this case Xz should drop as the compound 
is warmed up far above Ecp/kB since Az{c[,(jj = 0) be- 
comes negligible. This is inconsistent with the experi- 
mental finding for the garnet Yb3Ga50i2^*'. To under- 
stand the origin of this problem, we need first to spec- 
ify the Yb'^+ level scheme. Because Yb^+ is a Kramers 
ion, the states are at least doubly degenerate. The mean 
distance between the doublet ground state and the three 
closely spaced excited doublets is A//cb — 850 K— . Since 
the ground state is a doublet, Ecf = 0. In addition, 
(3 A ^ 1. Hence, from the RPA approximation we ex- 
pect an activated behaviour for Xz (T) outside the critical 
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regime (Tn = 0.05 K). But experimentally, from about 
0.4 up to 100 K, Xz is temperature independent. It dis- 
plays an activated behaviour only above '-^ 100 K. 

The breakdown of the RPA theory for Yb3Ga50i2 
could have been anticipated. In fact, the compound be- 
haves as a lattice of effective spins one-half coupled by 
the Heisenberg interaction up to ^ 100 K. This situation 
cannot be described by the RPA approximation we use 
because it assumes a non-zero single ion susceptibilitj«^. 
Only, if the temperature is increased to be such that the 
thermal energy is an appreciable fraction of the crystal 
field excitation energy, does Tier starts to matter and 
the material given in Sees. IIIll and IIVI becomes relevant. 
The Orbach relaxation then becomes visible. 

To model the quasi-elastic response and Xz for 
Yb3Ga50i2 at intermediate temperature, a finite expres- 
sion of the self-correlation function Kz{uj) for an isotropic 
magnetic material in the limit a; is needed. That 
function is given by the sum over the Brillouin zone of 

(q, Lo) . According to the fluctuation-dissipation the- 
orem, Kz{(\,u}) is proportional to the product of the 
static wavevector-dependent susceptibility, XzCq)? a-nd 
the spectral- weight function, i^z(q, oj). From the material 
gathered by Lovesey^, we could write a mean-field expres- 
sion for Xz(q) and ^^(q, w) in the Gaussian approxima- 
tion. Defining r^^oxc such that Az(w = 0) = 2( J^)/r2^cxc, 
we could identify Fz^exc by setting Kziuj = 0) equal to the 
sum over the Brillouin zone of Az(q, w = 0). However, 
the formulae for Xz(q) ^^^d -Fz(q, to) are approximate and 
the sum can only be done numerically. Here we propose 
a much simpler method to estimate Fz^oxc. We suppose 
the spectral function of the self-correlation to be Gaus- 
sian and use the expression of the second moment for 
paramagnets given in Ref. We get: 

r..cxc = V • (61) 

Interestingly, F^ oxc depends on the two- ion exchange in- 
teraction J (see Eq.HHl J — J[ij)) that we have hmited 
to the z nearest neighbour ions. Therefore it reflects the 
environment of the ions under study. F^ oxc is tempera- 
ture independent in the high-temperature limit we take 
and for the conventional magnets we consider, i.e. we 
do not discuss complicated compounds such as the frus- 
trated ones. 

The inclusion of the phonon-induced relaxation is eas- 
ily done. Referring to Eq. [55]for the definition of F^, 

Tz = Fz,cxc + Tz^ph. (62) 

In the same way, for the electronic relaxation we have 

Tz ~ rz,oxc + Tz.oi. (63) 

Fz,ei and Fz.ph are obviously obtained from the study of 
rizz(w). Eq. [nHhas been proven to provide a good fit of 
Xz{T) observed for YbaGasda^, but only for T > 0.4 
K. The increase of Xz for Tc < T < 0.4 K arises from 



the intersite correlations. The RPA provides a model for 
these correlations, but outside of the critical regime. 

In conclusion, the analysis of the neutron quasi-elastic 
linewidth and /iSR spin-lattice relaxation rate should be 
done, at least in a first approximation, with Eq. [54] or 
Eq. [55] and using for the single ion linewidth the expres- 
sion given at Eq. [S^] or Eq. [S3] 



VII. DISCUSSION AND CONCLUSION 

We have studied the symmetrised correlation functions 
as measured by neutron scattering and the /iSR spin- 
lattice relaxation rate. We have focused our work on 
the paramagnetic state of magnetic materials made of a 
regular lattice of rare-earth ions. The relaxation of the 
crystal-field levels has been assumed to be driven by the 
conduction electrons or the lattice vibrations. Mathe- 
matically, we have used the iterative Laplace transform 
method due to P. M. Richards. The interaction between 
the ions has been modelled with the random-phase ap- 
proximation. A phenomelogical modification of this ap- 
proximation has been proposed to cure its breakdown 
which comes into light when the crystal-field levels are 
not closely packed. Thanks to the point symmetry at 
the rare-earth site, for many compounds of interest, we 
have shown that at most two symmetrised correlation 
functions are relevant. 

Through the study of simple crystal-field models, the 
Orbach relaxation mechanism has been uncovered for the 
linewidth and Xz- We have argued for a single ion that 
the linewidth should vanish at low temperature if the 
relaxation is driven by the lattice vibrations. However, 
physically this limit cannot be reached because of the on- 
site fluctuations described by Fz.cxc in Sec. lVIl which set 
a lower bound on the linewidth and therefore an upper 
bound on Xz- 

As clearly seen from its derivation, the Richard's sys- 
tem of linear equations is obtained if the Hamiltonian de- 
scribing the interaction of the lattice field with the crystal 
field can be factorised into these two fields. This sets lim- 
its on the relaxation mechanisms which can be handled. 
For example, the Raman mechanism cannot be included. 
It involves two virtual phonons and it is described by 
treating the interaction Hamiltonian in second order—. 
It leads to a linewidth Fph c>c T"" with a large value for 
the exponent n, typically rt = 5 or n = 7. According to 
Electron Paramagnetic Resonance data, this mechanism 
can appear only at low temperature^. 

The modulation of the crystal field has been expanded 
in powers of the strain and we have kept the linear term. 
We could have gone a step further, keeping the term 
quadratic in strain. This would have introduced at least 
one new free parameter. Unless really justified for a par- 
ticular physical problem, we do not believe it is worth- 
while. 

Although our interest in this paper has been on the 
paramagnetic phase of a magnetic material, the method 
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we have discussed can be extended to its ordered mag- 
netic state. One has to account for the molecular field 
acting on the rare-earth ion by adding a Zeeman term to 
TicF as done in Ref. [27| . The RPA models the magnetic 
fluctuations due to the two-ion interaction. 

In conclusion, we have proposed a framework to anal- 
yse the linewidth as measured by neutron scattering and 
the /xSR relaxation rate for compounds with rare-earth 



sublattices. We take into account the electronic and 
phonon-induced relaxation of the crystal-field levels. We 
have proposed a phenomelogical modification of the re- 
sult of random-phase approximation which appears for 
some particular crystal-field level schemes. 

I thank Pierre Dalmas de Reotier for many discussions 
and his checking of equations. 
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